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We discuss the electromagnetic corrections to the ttN scattering lengths 
generated by minimal e. m. coupling from a knowledge of the low energy 
i/-) . expansion of the ttN elastic scattering amplitude as well as from the nucleon 

and A pole terms, all taken for purely strong interactions. We assume the 
heavy baryon limit; the e. m. and axial form factors and masses are assumed to 
have their empirical values, such that there is no free parameter. The different 
terms have a clear physical and intuitive origin. In particular, a large isospin 
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breaking contribution to the isoscalar term appears in the elastic charged- 
pion scattering lengths. We attempt a comparison to the results from chiral 



<]} [ effective field theory (EFT) with a physical interpretation of the empirical 

constants in that approach. The results are applied to the energy shift and 
width of the ir~p atom. 

1 INTRODUCTION 

In the limit of purely strong interactions the irN scattering amplitudes at threshold are 
fundamental quantities which enter into the discussion of various problems. They provide, 
for example, a basic test of the Tomozawa- Weinberg chiral relation for the isoscalar and 
isovector scattering lengths in the limit m n = 

a' = u/(8irF%) ~ 0.089 m" 1 ; a + = , (1) 

where F n = 93 MeV is the pion decay constant PI |2| and u = m n at threshold. The 
empirical isovector scattering length is the main ingredient and uncertainty in the forward 
dispersion relation, by which the irNN coupling constant 14.11±0.05 is determined^ etc.. 
The major precision source for these quantities are the remarkable measurements of the 
Is level shifts in pionic hydrogen and deuterium as well as the corresponding widths. In 
the case of pionic hydrogen, present experiments have, or will shortly achieve, a precision 
of 0.2% for the shift and 1% for the width |U|3] and these quantities convert in principle 
to similar precision for the scattering amplitudes. 
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There are 2 major ways by which one can approach the problem of determining the 
ir~p scattering length from pionic hydrogen data. The first is chiral effective field theory 
(EFT) [HI IE] Here one starts from an effective Lagrangian in the chiral limit and 
makes first a systematic expansion in orders 0(p n ) of momentum corrections. The e. m. 
contribution to order a and the strong symmetry breaking are then expanded as small 
additional corrections. This approach introduces empirical constants to summarize short 
range contributions and these must be determined from experiments. 

The other major approach, which I will follow below, is complementary and less ambi- 
tious |HlCn|. It starts from the empirical low energy expansion of the strong interaction 
amplitude in terms of energy and momentum. It is rather well explored including its prin- 
cipal dynamic features, although the experimental value of the scattering length needs 
further tuning. This system has well defined e. m. and axial form factors of the pion and 
nucleon (A) which are experimentally known as are the physical masses. Starting from 
this knowledge, we determine the e. m. corrections. The main guiding tool is minimal 
e. m. coupling or, in other words, current conservation. In addition, we use the same 
minimal coupling principle to determine the dispersive contributions from the radiative 
capture processes dominated by the nucleon and A isobar intermediate states. For this 
last problem, it is a considerable simplification to work in the heavy baryon limit and we 
will do so below. This approximation is also used in the EFT approach. 

The objective of this second approach is to obtain an intuitive picture of the corrections 
and the physical mechanisms of isospin breaking. It will become clear that kinematic 
considerations in the wide sense is the key to several of the effects. 

The plan is the following. I first discuss the nature and physics of the corrections 
to the leading order Deser-Trueman formula fTJ ^] , which relates the energy shift to 
the scattering length. I then discuss 'inner' e. m. corrections to the s wave scattering of 
charged pions from the nucleon at threshold and show that the dominant term is intimately 
connected to the well established p-wave irN scattering physics and the A isobar. This 
contribution can summarized by the EFT chiral parameter fi [HI EJj. Finally, I make 
a tentative comparison with the results of the EFT approach. 

1.1 Step 1: How to get scattering lengths from atomic energy 
shifts: removal of the external Coulomb field 

The s-wave threshold amplitude for the strong interaction has in the case of a single 
channel the low energy expansion 



where q is the momentum and bh are the scattering length and the range term re- 
spectively. Such an expansion does not by itself assume isospin invariance, but it means 
implicitly a short ranged interaction. 

The atom is highly non-relativistic. To leading order the strong interaction shift is 
obtained using the Is wave function for the point Coulomb field 
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4>Bohr(r) = 4>Bohr(0) exp(-amr) ~ </W r (0)(l - amr + ..), 
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where m is the reduced mass and amr << 1 over the strong interaction region. If 
we assume that the amplitude of Eq. (j2J results from a short ranged pseudo-potential 
('effective interaction Lagrangian'), the leading order Is strong interaction energy shift is 
simply the Born approximation term: 

Air 

eu = -—m 2 a h . (4) 
2m 

This is the time-honored standard reference shift referred to as the Deser-Trueman for- 
mula IHJII2] an d there are of course corresponding expressions for any ns state. We now 
generalize this expression. 

-Assume that the Coulomb field without strong interactions to be generated by the 
extended charge distributions of the pion and the proton. The corresponding e. m. form 
factor F em (q 2 ) = F n (q 2 )F p (q 2 ) is experimentally known. Its origin is in principle irrelevant, 
although the issue of a Dirac versus a Pauli form factor of the proton must be faced later, 
when we compare to results of EFT in the heavy baryon limit. At this stage the Pauli 
form factor is the relevant one.. 

-This problem is highly non-relativistic and well described using wave functions. Out- 
side the extended charge, the exact solution is the Whittaker function. 

-In the limit that the strong interaction (J2J) has a negligible range compared to that 
of the charge distribution, the corresponding strong interaction energy shift has an exact 
solution 0. 

-For the energy dependence of the scattering length we use the 'minimal e. m. coupling 
principle" or gauge invariance and replace the energy of a charged particle ui — > u— eVc(r). 
Another way of saying the same thing is that one should locally at the interaction point 
use the correct interaction energy in the Coulomb potential. 

There are now 3 separate physical effects jHj- First, the starting wave function at 
the origin is improved using the Coulomb potential corresponding to the joint charged 
distribution of the pion and proton, which is simpler to handle than the singular point 
Coulomb interaction. Second, the interaction does not correspond to the free scattering 
at threshold of Eq.(J2J), but corresponds to the energy shifted to that of the Coulomb 
potential at the interaction Third, to second order in the scattering length, the binding 
gives characteristic cusp term which is nearly independent of assumptions. Finally, we 
always include tacitly a small, model independent vacuum polarization correction, which 
is conceptually irrelevant in the present context. 

The extended charge changes the Coulomb wave function, chopping off the singular 
linear behavior of the wave function near the origin, which comes from the 1/r behavior 
of the Coulomb interaction. The wave function then varies normally with r 2 at the origin. 
The wave function at the origin is then 

<t>Bohr{0) -> <pBohM [1 - am(r) em + 0(a 2 )} (5) 

Here (r) em is the expectation value over the extended charge density. The corresponding 
change in the energy shift is 

4-7T 

Se ls;wf = -—(j)Bohr(0) 2 (-2ctm(r) em ) a h ~ (-0.9%) e h (6) 
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Table 1: Coulomb corrections in percent. The vacuum polarization contribution is in- 



cluded in the total correction. (From Ref.jSj) 





Extended charge 


Renormalization Gauge term 


Total 




-0.853(8) 


0.701(4) -0.95(29) 


-0.62(29) 


Sr 


-0.427(4) 


0.701(4) 0.50(23) 


1.02(23) 


S-K+p^-K + p 


0.853(8) 


0.72(5) -1.71(29) 


0.35(29) 



Detailed numerical values are given in Table 1. The wave function modification is not 
specific to the atomic bound state problem. It has a near identical counterpart in the 
elastic n + p scattering close to threshold, but the interaction is now repulsive. For a 
neutron target such Coulomb corrections are of course absent, we should change the sign 
of the correction term Eq. (|5|). but there is no change for a neutron target, of course. 

The energy change in the scattering amplitude (j2J depends only on minimal coupling, 
such that it happens for the scattering near threshold of a positive pion as well, but with 
opposite sign. 

Lu^uj-e(l + T 3 )t 3 V c (r). (7) 

For a short ranged limit we must use the potential energy at the origin. 

The range term in Ref. (J2J) corresponds to the experimental isoscalar and isovector 
range terms b + ± b~ for 7i T p interactions with the corresponding n~p energy shift: 

Air 1 
5e lS]gauge = -— <PBohr(0) 2 (-a(-) em )bl P ~ (-1.0%) e h . (8) 

We note first, that it does not matter, to leading order, whether this range term is 
derived in Eq. (J2J from terms of the momentum q or from the energy uo: the result will 
be identical [Hj, for the interaction is nearly on the mass shell. Second, the range term bh 
has its own physics and should not be taken to be proportional to the scattering length. 

The third contribution is the cusp effect generated by rescattering to second order in 
the scattering length by the binding in the Coulomb potential. It depends only weakly on 
the form factors. Its physics is that the incident wave in the interaction is renormalized 
by the scattering (effective field effect) as is known in the present context since a long 
time. It is accurately 

Seis;renorm = -8iraa 2 h [2 - 7 + log 2 a - (log mr) em ](j) Bohr . (0) 2 ~ (+0.7%) e h (9) 

These corrections are exact to order 0(a 2 ) in the limit of a short ranged strong interaction. 

A crucial term is the second one corresponding to the energy shifted amplitude depend- 
ing on the tt~p range parameter bh- One might think that the main contribution would 
follow from the uj dependence of the dominant isovector Tomozawa- Weinberg scattering 
length (jl]). which indeed appears to generate corresponding terms to next-to- leading or- 
der in the EFT expansion jHlElEl- This is not the case. This isovector range term is 
largely canceled by the nucleon pole term and the net effect is quite small. Instead the 
contribution is mainly generated by the isoscalar range term proportional to u 2 . This 
term is proportional to uj. 
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7T = 7T 



7T = 7T~ 



71 = 7T ± 7T = 7T ± 




p n N = p, n 7 N = p,n N = p,n 7 JM — p,n 

(a) (6) 

Figure 1: (a) the Kroll-Ruderman radiative capture process (b) diagrams of 7riV — > X7 — > 
ttN reactions generating isospin breaking in low-energy elastic irN scattering. 



It is a fair question to ask how accurate the expression for this energy correction 
is in practice. This can be inferred in several ways. An easy and intuitive estimate 
is obtained modifying the level displacement expression of Eq. (Jl| taking the interac- 
tion density to have a form factor F str (r) instead of a point interaction. This means 
that the Coulomb interaction should be averaged over the interaction region folding in 
this additional form factor. The corresponding interaction shift becomes Sei s;gauge cx 
/ Fst^q^F^q^Gp^q 2 ) / 'q 2 d 3 q. Using standard monopole and dipole form factors for the 
pion and proton with the scale of the p-meson and observing that the form factor of 
the Tomozawa- Weinberg term is also a monopole associated with the p-meson, indicates 
a typical small uncertainty from this source of about 0.15%, which is beyond present 
experimental uncertainties. 

Up to this point we have not used the heavy baryon approximation. This means that 
for the proton charge form factor we should here use the Pauli form factor, which includes 
the 1/Mp terms from the magnetic charge density and not the Dirac form factor. More 
about that later. 

All of these considerations are of course quite general and apply with small modifica- 
tions for any hadronic atom. They are easily generalized to states ns of higher quantum 
number n as well as to the width. It is also easy to generalize the result to a repulsive 
Coulomb interaction such as 7r + p, but the results will then concern the scattering length 
at threshold. Finally they can be generalized to the situation of two coupled channels 
The problem of extracting the hadronic scattering length from the atomic ir~p energy shift 
is therefore solved to a precision of about 0.1% on the level of the Coulombic interaction. 



1.2 Step 2. Transverse photons; a dominant isospin breaking 
mechanism at threshold 

The previous Coulombic contributions correspond to longitudinal photons. The question 
is then of the importance of transverse photons. In fact, these generate a large isoscalar 
correction, i. e., a term which is the same for all the elastic charged pion amplitudes as 
Dr. A. N. Ivanov and myself have recently shown [TU] . 

A guide to the importance of such terms is the dominant contribution of the Kroll- 
Ruderman radiative capture process rx~p — > 772 at threshold [H^ (see Fig. la), which 
experimentally gives a Is width of 8% of the strong interaction shift |16j, a huge number. 
The term we consider is the corresponding dispersive shift (see Fig. lb). 

The matrix element for radiative capture can, for example, be derived from the Partially 
Conserved Axial Current (PCAC) relation using minimal e. m. coupling (or directly from 
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the nucleon pole term). 

= -mlF^ w (x); d^d^± ieA„, (10) 

where A^ is the e. m. 4-vector potential. This statement corresponds to saying that 
the we have electric dipole (El) radiation due to the discontinuity in the current or, in 
other words, a kind of transition radiation. We calculate the dispersive contribution in 
the heavy baryon approximation. This is particularly convenient, since in this limit and 
at threshold the radiation comes only from the vertex itself in the Coulomb gauge and 
not from the pion and nucleon. 
The characteristic features are: 

-the transition is an axial one, such that its strength is well defined, 
-it is natural to use the axial form factor Fa{q 2 ), which is empirically well approximated 
by a dipole shape (see, e. g., [T7| ITHj) 

F A (q 2 ) = (1 + q 2 /Ml)' 2 with M A = (960 ± 30) MeV (11) 

- typical energy denominators (pim,,-)^ 1 appear from intermediate states with the sign 
switching due to crossing. 

-More precisely, the n~p amplitude contribution gives a large nucleon isoscalar term of 
3% in the limit m n = 0. The contribution is for the n~p case 

(i i m *\xf(ni)_ 3a 9Ar> f°° d PP F l(P 2 ) f , s 

V + M N ) dhs ~8^Ff Jo p-mr-iQ- 1 ' 

However, from the SU3 symmetry point of view, the A isobar and the nucleon are basically 
identical, but for the NA mass splitting and weight factors. It would be unnatural 
to include the nucleon only. In the limit of no NA mass splitting, the inclusion of A 
intermediate states gives a multiplicative factor 25/9 as compared to the nucleon term. 
The increases the previous 3% to 9%, an enormous correction! 

Why is this so large? The reason is that the scale parameter is the axial mass M A and 
not the pion one m n which gives a factor 7 enhancement with respect to naive expectations. 
When the NA mass splitting of ~ 2m n is brought in, it cuts the 7A term by 50% to a 
total isoscalar correction of about +6%. This is still very large, but we expect relativistic 
kinematic factors to cut it additionally to about 4.8%. 

1.2.1 The case of m n ^ 

In this situation new characteristic terms appear of the type cm^ In {itl^/Ma) + dm^, 
which are generated both from the nucleon and the A intermediate states. These terms 
depend only weakly on the exact value of Ma- In the particular case of the nucleon 
intermediate state, the term proportional to m n In m n has the identical coefficient to the 
one found to third order in chiral EFT by Gasser et al. |S]. If we expand the amplitude 
gfim) j n Eq_ (Tl2*j) in terms of the small parameter x = = m^/M^ , we have in this 
case: 

(l + ^ W-) = ^4 \-Ma -mJ\n^ + ^ + 0(^)1 . (13) 
V M N J J 87r 2 FH32 A w \ M A 12 K M A J\ v ; 



Isospin Violation in Threshold 7rN Scattering 



Table 2: Contributions to the ttN scattering lengths from dispersive radiative capture 





with nucleon and A intermediate states 


in units of 10 3 m, 


1 (from Ref. 


BDI). 




= 0,u A = 


(3.0(1)jv 7 + 5.3(2) A7 )t| 




= 8.3(3) t% 






= 0,uj a ^ 


(3.0(l) JV7 + 2.4(l) A7 )tl 




= 5.4(2) t\ 






^ 0,cu A = 


(2.6(1)jv 7 + 4.6(l) A7 )t§ +(-0 


8 Nl + 0.7 A7 ) t 3 r 3 


= 7-2(2) t 2 3 - 


O.H 3 r 3 






(2.6(1)jv 7 + 2.5(l) A7 )t§ +(-0 


8at 7 + 0.3 A7 ) t 3 r 3 


= 5.1(2) 


0.5t 3 r 3 



When the A isobar is degenerate with the nucleon, this nucleon term strongly canceled 
by nearly a magnitude, such that the the dependence on the pion mass becomes negligible. 
However, when the NA mass splitting is introduced in accordance with observation, the A 
contributions are quenched such that the contribution from the nucleon term is partially 
restored. A small term in the pion mass of about 50% of the value for the nucleon only 
survives. Numerical values for the different cases are given in Table 2. 



1.3 Step 3. E. M. isospin violation and Chiral EFT 

In the last few years a scheme has been developed to calculate strong and e. m. isospin 
violation to leading order for nN scattering using field theory methods based on an 
effective chiral Lagrangian [HJ |Hl HSj The amplitudes are calculated in a systematic 
expansion in powers of momenta (EFT). I will not enter into the details of this expansion, 
but only sketch a tentative comparison of some specific points. It is important to realize 
that certain of the predictions of such an effective field theory are specific and outside 
our present approach, while, on the contrary, our approach here generates terms of higher 
order in the EFT description than those presently considered. In the heavy baryon limit, 
the e. m. isospin breaking in the ttN threshold amplitudes are related to the the e. m. 
mass of the nucleon and the np e. m. mass difference in the EFT beyond the purely 
kinematic effects [6 . To next-to-leading order these relations can be expressed in terms 
of 3 constants fi^s- In the case of n ± p elastic scattering, which I chose for illustration, 
this gives the following relations. 



M em = _ e 2 F 2 



n 



fl + f3 



M em = _ e 2pl 



fi + h + h 



<T P = -2Mh ± \h] (14) 



Following Ref. jH] , I omit the term generated by the physical mass difference between the 
charged and neutral pion, which is of no concern in the present context. 

It is of considerable interest to attempt to identify our results in the EFT expansion. 
Can we match our previous results to its explicit terms? A problem occurs in view of 
the basic difference with our approach. For example, in the discussion of the Coulombic 
terms, we start from the scattering amplitude and form factors as they would result in 
a description to all orders, but in the absence of e. m. interactions; in this sense we 
include physics not presently included in the EFT approach. We must therefore make 
sure we compare comparable quantities. First, we must use form factors in the heavy 
baryon limit. This means that we must use the Dirac form factor F p (q 2 ) and not the 
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Pauli form factor G p (q 2 ); the charge distribution generated by the magnetic moment 
of order M^ 1 should be omitted. We should also to this order omit the wave function 
correction am(r) em a n - p of Eq. (jSJ, for the mass scale is set by the p-meson and it is of 
order am/m^a^ or formally of 4th order in the EFT expansion 1 . The relevant term for 
the comparison is generated by Eq. (J5j). Here the dominant contribution comes from the 
isoscalar range term b + . It generates a contribution of order ma{r) em b Jr , which is of 4th 
order in the EFT expansion and outside the present EFT discussion; it will require an 
additional EFT constant. Similarly, the pole term in the isovector range term vanishes in 
the heavy baryon limit. 



1.3.1 The EFT constant / 2 , the np mass difference and the Coulomb interaction 

In the heavy nucleon limit, the constant f 2 in Eq. (pPfj) describes the np e. m. mass 
splitting, which then results from the Coulomb self energy of the proton with the Dirac 
form factor: (M p — M n ) em oc f 2 oc ~ J d 3 q[F p (q 2 )} 2 / q 2 . The corresponding effect in the 
7r c p scattering must be taken in the same limit with the Dirac factor. To leading order, 
the strong scattering amplitude is given by the Tomozawa- Weinberg term (TjQ), which is 
linear in u. The minimal coupling procedure generates an isospin-breaking contribution: 



n em - 1 n I r )t 2 BVc ^ MV 
a*± P ~ 2^ i + r 3^3 ( 87rF 2) ' (l0) 



where eVc(0) oc a J F n (q 2 )F 2 (q 2 )/q 2 d 3 q. Compared to the proton e. m. self energy f 2 
above, the similarity is striking and numerically the expressions are equal to about 15%. 
However, if the shape of the strong interaction is included with an interaction form factor, 
one obtains equality if 

F str (q 2 )F«(q 2 )F p (q 2 ) ~ F 2 (q 2 ) (16) 
or, since both the nucleon and pion form factors are governed by the p-meson mass 

F str {q 2 )~F«{q 2 )~(l + q 2 /m 2 p )-\ 

We here obtain unexpectedly a quantitative result for the isovector interaction form factor. 
The Tomozawa- Weinberg term is frequently associated with a monopole interaction with 
a range given by the p meson mass (1 + q 2 /m 2 )~ 1 j20], while both the nucleon and the 
pion have form factors closely approximated by a dipole, respectively a monopole, form 
factor with the p-meson mass. 

This is suggestive, since it indicates that, indeed, the EFT f 2 coefficient comes from the 
Coulomb interaction. On the other hand, the EFT does not presently give the shape of 
the form factor. The last relation implies however that there is an intrinsic link between 
the proton, pion and the Tomozawa- Weinberg interaction form factors in EFT. 

A proper account for the Coulomb interaction requires realistically that one goes beyond 
the present level of the EFT expansion and accounts for the Pauli form factor Ge- p of the 
proton as well as for the terms generated by the isoscalar range term. 



1 Since its coefficient is large, its order is unclear; its magnitude corresponds to 3rd order. 
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1.3.2 The EFT constant fi and the axial form factor 

The dispersive contribution from intermediate 7-/V(A) states to the scattering length as 
given in Eq. (Tl2"j) is isoscalar in the charged pion sector in the limit of a vanishing pion 
mass. This is exactly the symmetry property of the contribution in EFT by the next-to- 
leading order constant f\. This constant also appears as a part of the e. m. neutron mass 
M^ m , but then it always comes in the combination j\ + f% (see EqfJJ])). These terms 
cannot be physically separated [E]. Here we have an interesting situation. The neutron 
e. m. mass is expected to be quite small compared to the np e. m. mass splitting of 
about 0.8 MeV and to be dominated by the magnetic self energy. It then vanishes in the 
heavy baryon limit with a likely uncertainty of ±(0.1 ~0.2) MeV This corresponds 

to a value F%\f\\ ~ 1 MeV, while our value with the physical A isobar included gives 
F^fi = —26(1) MeV, which is over a magnitude larger. Dimensional estimates inside 
of EFT give intermediate estimates F%\fi\ = 6 MeV and 12 MeV [SJIE]- There appears 
therefore numerically to be little relation of this parameter with the neutron e. m. mass, 
which suggests a massive cancellation between the EFT constants /i and ^3. Such a 
cancellation occurs explicitly in a model evaluation inside a heavy quark model [22], 
which has certain similarities with our description. The connection of the constant j\ to 
the nucleon e. m. mass is therefore tenuous and of little practical importance. 

2 Conclusion 

Isospin violation in ttN elastic scattering has been shown here to have well determined 
contributions originating in the Coulomb field of the extended charge with little model 
dependence. These corrections are general and involve terms beyond present EFT ap- 
proaches. To leading order in the strong interactions corresponding terms have coun- 
terparts in EFT in the heavy baryon approximation. In addition, terms with r yN(A) 
intermediate states give rise to important isoscalar isospin breaking terms which can be 
looked at as model descriptions of the unknown EFT constant j\. For a non- vanishing 
pion mass, the same mechanism generates small, model insensitive, isospin breaking in 
the isovector interaction dependent on the pion mass. 

An important finding is that the isospin breaking is small in the isovector amplitude. 
This is consistent with the finding of Meissner et al. jSj, but it is in violent disagreement 
with the important violation reported by Matsinos [2*3] , which depends only on isospin 
breaking in the isovector amplitude, although the author does not explicitly state so. 
The Matsinos results have also been shown to be grossly at variance with the empirical 
scattering lengths deduced from pionic hydrogen and deuterium [3] . 

Although therefore the chiral EFT constants / 1)2 can be rather well understood from 
physical effects, this is not the whole story. Our terms include the main part, but not 
all, of the isospin breaking mechanisms. In addition to the effects discussed here, the 
effective field theory generates generic isospin breaking beyond our considerations. These 
additional terms appear in charge exchange and n°N elastic scattering. They are are 
conceptually important non-trivial chiral predictions. Numerically, such contributions 
are of the order of 1% and of similar magnitude as the terms we consider. Obviously 
a quantitative test of the predictions of EFT requires that that these intrinsic terms be 
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reliably separated from those which are generated by the general mechanisms discussed 
here. A meaningful analysis will require both approaches. 
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